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ABSTRACT 


A FORTRAN COntrol Program for Engineering Synthesis 
(COPES) has been developed for solving engineering design 
problems. The program maximizes or minimizes a numerically 
defined objective function subject to a set of inequality 
eoustraint.functions. COPES ‘uses the optimization program, 
CONMIN, which includes the conjugate direction method of 
Fletcher and Reeves for unconstrained function minimization 
and a modification of Zoutendijk's method of feasible direc- 
tions for constrained function minimization. Additionally, 
approximation techniques are available for use in optimiza- 
tion, and trade-off studies may be performed. A simple 
design example demonstrates the program capabilities. 
Programming guidelines are presented followed by sample 


input data and output for each program option. 


I. INTRODUCTION 


Most design processes require the minimization or maximi- 
zation of some parameter which may be called the design ob- 
jective. For the:design,to.be ;acceptable, .it_must satisfy a 
variety of physical, aesthetic, economic and, on occasion, 
political limitations which are referred to here as design 
constraints. While part of the design problem may not be 
easily quantified, most of the design criteria can be des- 
cribed in numerical terms. 

To the extend that the problem.can be stated in numerical 
terms, a computer program can be written to perform the 
necessary calculations. For this. reason,. computer) analysus 
is commonplace in most engineering organizations. For example, 
in structural design, the: configuration,, materials. and loads 
may be defined and a finite element analysis computer code 
is used to calculate stresses, deflections and other response 
quantities of interest. If any of these parameters are not 
within the prescribed bounds, the engineer may change the 
structural member sizes and re-run the program. The computer 
code therefore provides only the analysis of a proposed design, 
with the engineer making the actual design decisions. This 
approach es design, which may be called computer-aided design, 


is commonly used today. 


Another common use of analysis codes is in trade-off 
studies. For example, an aircraft trajectory analysis code 
may be run repetitively for several payloads, calculating the 
aircraft range, to determine the range-payload sensitivity. 

A logical extension to computer-aided designe is- fully 
automated design, where the computer makes the actual design 
decisions, or to perform trade-off studies with a minimum of 
man-machine interaction. The purpose of the COPES program 
is to provide this automated design and trade-off capability. 
The user must provide a FORTRAN analysis program for analysis 
of the particular problem being considered. This analysis 
program is written according to a simple set of guidelines 
so that it can be easily coupled to the COPES program for 
automated design synthesis. 

This document describes the capabilities of the COPES 
program and its usage. A simple design example is first pre- 
sented to demonstrate the program capabilities. Guidelines 
are given for writing analysis codes which can be coupled 
directly to COPES. Finally the data organization is outlined 
and sample data is presented. 

The most commonly used option of the COPES program will 
be for design optimization. Two approaches are available 
for this purpose. The first (NCALC = 2) is direct optimiza- 
tion of the function by the CONMIN optimization sub-program 
[Ref. 1]. An alternative to this is through the use of 


approximation techniques (NCALC = 6, Ref. 2). This second 


Option, 1S. usually, more, effivetent. for problems of underas 
design variables, but.which requires costly. analysis, espes 
cially, when, multaple optimizations are \to be performed: 
DESIGN EXAMPLE 


P=10, ie ib: 


Assume it,.is required to 
design the cantilevered beam 


shown, in Fig. J... The jobjgec- 
tive as to) find che minimum 


volume of material which will 


Support the concentrated load. |— L-200 in. — | 
Thavt, 1S F=30x10° psi 
Minimize Volume = BeHeL (a) 


Figure, d —. CantilLevyened 
Beam 


The bending stress in the beam must not exceed 20,000 psa 


_ Mc _ 6PL 
= 7~ = 7H < 20,000 (2) 


The shear stress must not exceed 10,0007 ps1; 
Vis ee peo tele O00 C3) 


and the deflection under the load must not exceed one inch; 


Spi 2 Ape 
§ = eet See ec 1.0 (4) 


Additionally, geometric limits are imposed on the beam size 


so that; 
0:.5, < B. <, 5.0 (5) 
Ie 0 sare 20°20 (6) 
H/Bg< 140'..0 CD 


Equation (1), the design objective, provides a measure of 
the efficiency of the design, while Eqs. (2-7) define the 
criteria which the structure must satisfy, and are referred 
Oo: as constraints. 

This design problem may be stated in general form; 


Minimize F(X) 


Subject to: 

G(X) < 0 ab, = am 
x* ae eee 

i-°i- ‘i 


where X is a-vector containing the design variables B and H 
and G(X) are the constraints defined by Equations (2-7). 
There are eight constraints on the design. The objective 
and constraints are functions of the design variables. 

While this problem is straightforward, its solution is 
not trivial because it is not known which constraints will 
be critical (i.e. which G(X) = 0) at the optimun. 

One design approach is to assume that the bending and 
the displacement constraints are critical, so that Equations 
(2) and (4) are equalities, and then solve the two simulta- 
neous equations for the design variables B and H. The 
remaining constraints are then checked to be sure they are 
satisfied, and if not, they are used to obtain a new design 
satisfying all constraints. When a design is obtained where 
two G5 (X) = 0 and the remaining G(X) < 0, PE. 1S am accept 


able design. 


The design obtained here yields a volume of 6750 and 
satisfies all constraints. Note, however, that the objec- 
tive function of Equation (1) played’ no part in” the ‘chotee 
of design variables so that, using this approach, there is 
no assurance that a minimum volume design will be produced 
[Ref. 3]. More importantly, it is desirable to devise 
techniques for optimum design of systems which may be de- 
fined by more than two variables and by much more complex 
analysis, where the optimum design cannot be determined by 
inspection. 

The, COPES. program provades this general capability. by. 
the, use of, the optimization, program CONMIN [Ref.. 1]. Tosuse 
this design capability, a FORTRAN code must be provided 
which, will calculate. the various parameters: In writing the 
analysis, code; 1) it is written in, subroutine form with 
SUBROUTINE ANALIZ (ICALC) as the main routine, 2) it is seg- 
mented into INPUT, EXECUTION AND OUTPUT and 3) all parameters 
which may be design variables, object functions or con- 
straints are contained in a single labeled common block called 
GLOBCM. 

To. demonstrate the simplicity with. which a desien- 
oriented analysis code can be written, the following FORTRAN 
subroutine was produced for the analysis of the cantilevered 


beam in Figure |. 
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SUBROUTINE ANALIZ (ICALC) 
COMMON /GLOBCM/ B,H,VOL,BSTRES,SHRSTR, DELTA ,HB ,E,AL 
ie (LCALCG.GT.1) GO-TO 16 
C --- INPUT OR INITIALIZATION. 
B=2.5 
H=10. 
P=10,000. 
E=30.E+6 
AL=200. 
WRITE (6,30)AL,P,E,B,H 
RETURN 
@=--, EXECUTION. 
10 CONTINUE 
VOL=AL*B*H 
BSTRES=6 .*P*AL/(B*H*H) 
SHRSTR=1.5*P/(B*H) 
DELTA=4.*P*(AL**3) /(E*B* (H**3) ) 
HB=H/B 
IF (ICALC.LT.3) RETURN 
C --- PRINT RESULTS. 
20 WRITE- (6,30) AL,P,E,B,H 
WRITE (6,40) VOL,BSTRES,SHRSTR, DELTA ,HB 
30 FORMAT(/////5X, 17HCANTILEVERED BEAM//5X, 8HAL =,F9.3/5X, 
* 8HP =,E12.5/5X, 8HE =,E12.5//5X, 8HB = ,F9.3/5X, 
* 8HH =,F9.3) 
40 FORMAT(/5X,8HVOL  =,F9.3//5X,8HBSTRES=,E12.5/5X, 
8HSHRSTR= 
* £12.5/5X, 8HDELTA=,E12.5/5X, 8HH/B =,F9.3) 
RETURN 
END 


This routine may be executed as a simple analysis program 
by the following main program; 
C MAIN PROGRAM TO EXECUTE SUBROUTINE ANALIZ. 
DO 10 I = 1,3 
10 CALL ANALIZ (1) 
STOP 
END 


Moreover, ANALIZ can be coupled directly to the COPES 
program to perform this same function, or to perform optimiza- 


tion or trade-off studies. 


This subroutine was coupled to COPES and an optimization 
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was performed (NCALC = 2) to yield the following design; 


CANTILEVERED BEAM 


L = 200.00 

P = .10000E+05 
E = .30000E+08 
B v5 te) 

H = 18.16 

VOL = 6607.61 
BSTRES =. 2.000:2E+0'5 
SHRSTR = .4540ZE+03 
DELTA =. 9790 60E+00 
H/B ="9.9'8 


The parameters L,-P,"E, Boand.H were’ input? tof the pre- 
gram and the remaining parameters were calculated. The design 
variables B and H were changed during the optimization process 
to. obtain this: design’. 

This design was achieved with 45 calls to ANALIZ with 
ICALC = 2 (45 analyses). This design could surely have been 
found with fewer analyses were it performed bythand calcula-— 
tions. However, once having written the analysis subroutine, 
numerous other designs may be obtained for different materials, 
loading, or design stresses with only minimal effort. 
Furthermore, the design obtained here of vollume!= 6607.6 1s 
very near the theoretical optimum of 6603.9, while that whven 
Satisfied stress and displacement constraints simultaneoushy 
was not. 

This very simple design example underscores the power of 
numerical optimization techniques and the ease with which they 
may be applied. The key to efficient use of the COPES program 


is the requirement that the ANALIZ code be written in a standard 
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format so that it may be coupled to COPES without modifica- 
tion. The following section contains guidelines for writing 


a design-oriented analysis code. 


id 


II. PROGRAMMING GUIDELINES 


In developing any computer code for engineering analysis, 


Zt. /1S\ prudent, to wratesthe code in such a way that a4 1s 


easily coupled to a general synthesis program such as COPES. 


Therefore, a general programming practice is outlined here 


which in no way inhibits the use of the computer program in 


its traditional role as an analytic tool, but allows for same 


ple adaptation to COPES. This approach is considered good 


programming practice and provides considerable flexibility 


of design options. Only five basic: rules must be followed: 


iP. 


i 


Write the code in subroutine form with the primary 
routine called as; SUBROUTINE ANALIZ (ICALC) The 
name ANALIZ is compatable with the COPES program and 
ICALC is a calculation control. Note that subroutine 
ANALIZ may call numerous other subroutines as re- 
quired to perform the necessary calculations. 

Segment the program into INPUT, EXECUTION and OUTPUT. 
The calculation control, ICALC, will determine the 


portion of the analysis ‘code to be executed. 


ICALC = 1; the program reads all data required to perform 


the analysis. Also, any initialization of constants which 


will be used repetitively during execution is done here. 


This initial input information as printed here tor, later =sne- 


ference and for program debugging. 
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ICALC =.2: the program performs the execution phase of 
the analysis task. No data reading or printing is done here, 
except on user-defined scratch disc. Data may be printed 
here during program debugging, in which case it should be 
controlled by a print control parameter which is read during 
input. In this way, this print may be turned off after the 
program is debugged, but may be used again during future 
program expansion and debugging. The reason that printing 
is not normally allowed during execution is that when 
optimization is being done, the code will be called many 
times with ICALC = 2, resulting in voluminous print. 

ICALC = 3: the results of the analysis aré printed: 
Also the essential input parameters which may have been 
changed during optimization should be printed here for easy 
reference. In some cases, so much information is generated 
when ICALC = 2 that it would be inefficient to store it 
internally or on disc for printing when ICALC = 3. In this 
case, it may be desirable to actually execute when ICALC = 3 
with a print code turned on to print his intermediate infor- 
mation. However, this approach should be avoided because 


it requires an additional complete execution of the program. 


In summary, when; 


ICALC = 1 Read input data. 
ICALC = 2 Execute the analysis. 
ICALC = 3. Print the results. 


iS 


III. Store all parameters which may be design variables, 
abjectane functions or constraints in a single 
labeled common block called GLOBCM. The order in 
which itheytare ‘stored is arbitrary. A lrsting of 
the COPES program should be checked to see how many 
parameters my be stored in GLOBCM (the dimension of 
ARRAY)/2 > “Inatialedistributiron of CORES allows fom 
1500 parameters. 

IV. During execution or output, no parameters which ame 
read during input should be updated. For example, 
if variable X 1s anitaalazed during input; the 
execution segment must not update X such as X = X@ 
3.2. Instead “a new variable, Y = xX + 372° shouldbe 
defined. 

V. Write all programs in standard language, avoiding 
machine dependent capabilities such as multiple 
entry point (IBM), DEFINE statements (UNIVAC) and 
seven letter FORTRAN names (CDC). “While this guide= 
line is not essential to the use of the analysis 
code within the COPES program, it makes the analysaie 
code much more transportable between different 
computer systems, a capability which easily justi- 
fies a slight reduction in efficiency on a given 
machine. 

Adherence to these guidelines not only leads to a more 


readable and machine independent computer code, but allows 
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this code to be coupled to the COPES program without 
modification. 

Having written the analysis code, it may be executed 
either with a simple main program or within the COPES program 
to perform the analysis. To insure that guideline IV is 
followed, the following main test program is recommended. 
Note that this program calls ANNALIZ twice with ICALC = 2 and 
ICALC = 3, to show that the same result is obtained repetitively 


C MAIN PROGRAM TO CHECK SUBROUTINE ANALIZ, 
C READ, EXECUTE AND PRINT. 
DO: 10 TCALC2>..1,--3 
10 CALL ANALIZ (ICALC) 
EXECUTE AND PRINT AGAIN TO BE SURE THE RESULTS 
C DO NOT CHANGE. 
DO 20 ICALC = 2, 3 
20 CALL ANALIZ (ICALC) 
STOP 
END 


This program was executed with the ANALIZ subroutine for 
the cantilevered beam example to yield the following result: 


CANTILEVERED BEAM 


AL = 200.00 

P = .10000E+05 
E = .30000E+08 
B = 2.00 

H = 5.00 
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CANTILEVERED BEAM 


AL = 2010/00 
Ie =e 710000B705 
E = .50000E+08 
B =" 2.00 

=e oy (0)(0) 
VOL = 2000.00 
BSTRES = .24000E+06 
SHRSTR =. LSOOOETO4 
DELTA = .42667E+02 
H/B = 2.50 


CANTILEVERED BEAM 


AL = 200.00 

4 = .10000E+05 
E = .30000E+08 
B = 2.00 

H = 5.00 

VOL = 2000.00 
BSTRES = .Z4000E+06 
SHRSTR = .15000E+04 
DELTA = .42667E+02 
H/B =" (2.50 


This design was used as the initial design for the optimi- 
zation presented previously. Note that, while the volume here 
is lower than the optimum, the bending stress and displacement 
each exceed the imposed limits by more than an order of 
magnitude. 

Once the analysis code has been written, it can be coupled 
to the COPES: program without modification. ~lf it 1s desired 


to perform a simple analysis using  COPES, only three data 
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cards are required for the COPES program, namely a TITLE card, 
a control parameter, NCALC = 1, and an END card. If the 
optimization or parametric analysis (sensitivity) capabilities 
of COPES are to be used, additional data must be read. This 
data will identify which parameters in the global common block, 
GLOBCM, are used. To set up the COPES data, the user must 
have a basic understanding of how the data in the global com- 
mon block 1s accessed by COPES. This is outlined in the 


following section. 


aE 


III. DATA MANAGEMENT 


In order to perform design operations, the COPES program 
must access the data in common block GLOBCM. This is done by 
defining the location in GLOBCM where a specified parameter 
resides. For example, consider the common block for the can- 
tilevered beam; 

COMMON/GLOBCM/B,H,VOL,BSTRES ,SHRSTR, DELTA,HB,E,AL 

The volume of material, VOL, is the third parameter in 
the common block; that is, it resides in location 3, referred 
to as the global location number. Similarly the bending 
stress, BSTRES, is in global location 4 and the beam width as 
in global location 1. Thus, the parameters are referredito 
by their respective location numbers in global common. 

For convenience in preparing data for the COPES progran, 
a simple "CATALOG" of parameters may be defined. For the 


cantilevered beam, this catalog would be; 


GLOBAL FORTRAN 
LOCATION NAME DEFINITION 
i B Beam width 
Z H Beam height 
3 VOL Volume of Material 
4 BSTRES Maximum bending stress 
5 SHRSIR Maximum shear stress 
6 DELTA Deflection under the load 
7 HB Ratio, H/B 
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GLOBAL - FORTRAN 


LOCATION NAME DEFINITION 
E Young's modulus 
5) AL Length of beam 


As another example, consider a global common block con- 
taining arrays; 
COMMON/GLOBCM/A, Y(10), Q, C(2,2), H 


The variable catalog for this common block is; 


GLOBAL FORT RAN 
LOCATION NAME DEFINITION 
i A Area 
Z Y(10) Vector of y-coordinates 
LZ Q ° 
13 C(Z,2) . 
L7 H ete, 


The dimensions are given with the FORTRAN name as a 
reminder that the parameter is an array. In this case, the 
third parameter in the Y array is in global location 4. 
Remembering that arrays are stored column by column the 
C(1,2) array location is in global location l15. 

It will be seen that identifying parameters according 
to their location in GLOBCM provides a great deal of flexi- 
bility in using the COPES program for design. 

Appendix A contains blank forms for writing the variable 


catalog for the user's particular problem. 


ZL 


IV... COPES TERMINOLOGY 


The copes program currently provides six specific 


capabilities: 

ly “Simple “analysis, just as sif COPES@wasi not usede 

2. Optimization - Minimization or maximization of one 
calculated function with limits imposed on other 
functions using the CONMIN subprogram. 

3. Sensitivity analysis — the ettiect, of ‘changing jome 
or more design variables on one or more calculated 
functions. 

46 Two-variable functaon spaces analysis for vail 
specified combinations of two design variables. 

5. Optimum sensitivity - samevas sensitivity sanaljysus 
except at each step, the design optimized with 
respect to the remaining independent design vari- 
ables. 

6. Approximate optimization - optimization using 


approximation techniques. Usually more efficient 
than standard optimization for up to 6 design 
variables or if multiple optimizations are to be 


performed. 


In defining the data required ito execute jthe ,\COPES mro- 


gram, 


the following definitions are useful: 


Design Variables - Those parameters which the optimiza- 


tion program is allowed to change in order to improve the 


design. 


Design variables appear only on the right hand side 
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of equations in the analysis program. COPES considers two 
types of design variables, independent and dependent. If 
two or more variables are always required to have the same 
value or be in a constant ratio, one is the independent 
variable while the remaining are dependent variables. For 
example, if the height is required to be 10 times the width 
of the cantilever beam, B would be the independent variable 
while H would be the dependent variable. 

Objective Function - The parameter which is to be mini- 
mized or maximized during optimization. Also the parameters 
calculated as functions of specified design variables during 
a sensitivity or two-variable function space study. Objective 
functions always occur on the left side of equations, unless 
the objective function is also a design variable (the beam 
height may be minimized as an objective function if it is 
also a design variable. In this case the minimum height is 
found for which no constraints are violated). An objective 
function may be linear or non-linear, implicit or explicit, 
but must be a function of the design variables to be meaning- 
ful. 

Constraint - Any parameter which must not exceed speci- 
fied bounds for the design to be acceptable. Constraint 
functions, always appear on the left side of equations. Just 
as for objective functions, constraints may be linear or 


non-linear, implicit or ‘explicit, but must be functions of 


the design variables. 


25 


Constraint Set - A group of constraints which appear 
consecutively in the global common block and which all have 
the same limits imposed. “This is a convenience which allows 
several constraints to be identified with a minimum of data. 

Global Common - Common block GLOBCM containing design 
information. 

Global Location - Location of a particular parameter in 


GLOBCM. 
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V. COPES DATA 


The COPES program reads data from unit 5 and writes output 
on unit 6. Units 20 and 40 are used as scratch files. The 
scratch file numbers may be changed by changing two cards at 
the beginning of the COPES program. 

In order to execute the COPES program it is necessary to 
provide formatted data for COPES, followed by data for the 
ANALIZ program which is coupled to COPES. This section de- 
fines the data which is required by COPES. The data is seg- 
mented into "BLOCKS" for convenience. All formats are alpha- 
numeric for TITLE, and END cards, F10 for real data and 110 
for integer data. 

The COPES data begins with a TITLE card and ends with an 
END card. This is followed by data to be read by the user 
supplied subroutine ANALIZ or when ICALC = 1. 

Comment cards may be inserted anywhere in the COPES data 
stack prior to the END card, and are identified by a dollar 
sign ($) in column 1. 


Data coding forms are provided in Appendix B. 


VI. UNFORMATTED DATA INPUT 


While the user's sheet defines COPES data in formatted 
fields of ten, the data may actually be read in a simplified 
fashion by separating data by commas or one or more blanks. 

If more than one number is contained on an unformatted data 

card, a comma must appear somewhere on the card. If expo- 

nential numbers such as 2.5+10 are read on an unformatted card, 

there must be no embedded blanks. Unformatted cards may be 

intermingled with formatted cards. Real numbers on an un- 

formatted card should have a decimal point. 

Examples: 

Unformatted data; 

Di 5: Lapelivegd ulnes Oct 2iOus 0 at couaek 

Syl aoahts: Senta: Ota Oke n5) 2 

57 AS 1 0F20 5 eee 

5 Fan ls dig Age OFZ0s, 0. = 5 n2 

Equivalent formatted data; 

col> 10 20 30 40 50 60 70 80 
fi Led  L.0HZ20 0 -$.1 


Unformatted data; 


26 


Equivalent formatted data; 


col>.- 10 20 30 
Z 
Note: This data has been right 
7 3 justified. 
Z 35 Note: This data contains no commas, 


so it is assumed to be format- 
ted already. 


Unformatted data; 
25 5544 550,7,849, 10,11 


Equivalent formatted data; 


col> 10 20 30 40 50 60 70 80 

a OY 
i 2 5 4 5 6 7 8 
2) 10 EL 


Note that two formatted data cards are created here. 
Unformatted data; 

235545550 

7,55, 9,10,11 


Equivalent formatted data; 


col> 10 20 30 40 50 60 60 80 
° 
fi 8 9 10 171 


Note that the above two examples do not produce the same for- 


matted data cards. 
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VII. SAMPLE COMPUTER SOLUTIONS 


Sample solutions are presented here for the various COPES 
program options defined by the parameter NCALC in DATA BLOCK 
B. The ANALIZ routine given in the design example section was 
used to produce the results given here. Note that SUBROUTINE 
ANALYZ does not read data, so only COPES data is required here. 
In the usual case where data is read as input to the analysis 
routine, this data would follow directly after the COPES. 'END' 
card. 

| The output from a COPES program execution includes a title 
page followed by a copy of the input data. Then the required 
program executions are performed and final output information 
iS. printed. Figures 2-7 Contain output for all options of 


COPES as follows: 


NCALC FIGURE PROBLEM SOLVED 
iI Z Analysis only 
2 D Optimization by CONMIN 
3 4 Sensitivity study 
4 ) Two-variable function space study 
5 6 Optimum sensitivity 
6 7 Optimization using approximation 


techniques 


ao 


Note that data for the COPES program is read depending on 
the value of the parameters NDV, NSV, N2VAR and NXAPRX in DATA 
BLOCK B. The actual program execution is determined by the 
value of NCALC. ‘Therefore; ’ data may'’be read’for all*program 
options even though only one option of the program is to be 
executed. Figure 8 is a copy of the combined data for all the 
previous examples. This data may be used for any program op- 
tion simply by changing the’ value-of NCALG: -+Pigure/9Mismee 
copy of this same data using the simplified data input mode 


and Figure 10 1s the same data without the comment cards’ 
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